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A strength function method is adopted to describe a coupling between 
electric and magnetic modes of different multipolarity. The collective vibra- 
tions are analysed for a separable residual interaction in the framework of 
the random-phase approximation. The coupling between M2 and El giant 
resonances is considered as an illustrative example. 
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I. INTRODUCTION 

The probe of nuclear structure via electromagnetic interactions provides important in- 
formation on the nuclear wave function and dynamics of nucleon-nucleon interactions. The 
nuclear response upon different external fields is measured by means of electromagnetic 
transitions between different quantum states. The transition amplitude at the emission (ab- 
sorption) of a photon of a given multipolarity is proportional to matrix elements of multipole 
moments. Since in nuclear photo-processes the wave length A ~ 1/k of the photon is larger 
then the nucleus radius R (the long- wave approximation kR << 1) the multipole moments 
take the form |JTJ] 
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where is the nucleon effective charge of multipolarity A, gs^ft(i) and g\ X }ff(i) are the 

corresponding effective spin and orbital gyro magnetic ratios, respectively, and hn = ■^r is 
the nuclear magneton. In this approximation the total decay probability from the quantum 
state \i > to the quantum state \f > is determined as 
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where the reduced transition probability for the electromagnetic transition of type X 
(X = E or M) and multipolarity A is defined as 
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Here the B(EX,i — > /) is given in e 2 fm 2X , B(MX,i — > /) is given in /j, 2 N fm 2X ~ 2 , the 



transition energy E 1 is defined in Me V and the total decay probability T ir __>f(E 7 ) is given in 
1/s. 

At high excitation energy the density of nuclear states becomes so large that the descrip- 
tion in terms of transition probabilities between individual states loses its practical sense. 
Strength function techniques are known to represent a very valuable tool for studying the 
electromagnetic properties of nuclei in energy regions with high level density For 
example, defining the strength function as 

S(XX, £ 7 ) = J2 B(XX, gr - u)6{E y - u„), (4) 

V 

we can describe the distribution of strength of the electromagnetic excitation of the ground 
state over a range of the available excitation energy E 1 . Consequently, the strength function 
can be used for the analysis of the photo-absorption cross section cr(E^) 
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Here, S(EX,i — > f), S(MX,i — >• /) are given in e ^ v , N MeV — , respectively. The compari- 
son of the strength functions extracted from photo-absorption experiments with theoretical 



estimates can testify the basic model assumptions used for the calculations. Notice, that in 
the expression Eq. (|5|) there is a restriction due to only the parity selection rules. Therefore, 
such excitations contain contributions from different electric and magnetic operators which 
are able to create the excitation with a given parity. 

In axially deformed nuclei rotational bands are built on the intrinsic states characterized 
by the angular momentum projection K and the parity. In fact, such intrinsic states can be 
created by different electric and magnetic operators with all possible multipolarities in the 
nuclear Hamiltonian. However, there is a restriction on the parity quantum number again. 
For example, the state with K n = 1~ can be created by El, M2, E3, M4... operators. In 
standard approaches the coupling between different electric and magnetic modes allowed 
by the parity selection rules is neglected ]2]||, however, there is no justification for this 
assumption. Based on the strength function method ||, we propose the generalisation 
of this approach |4j] which allows consistently treat the contribution of different electro- 
magnetic modes to the strength of excitations (de -excitations), in particular, in region of 
giant resonances. 



II. THE MODEL 

For the analysis of collective and single-particle degrees freedom of nucleus we use the 
the phonon-l- rotor model described in details in |4j]. The model Hamiltonian is 

H H S p -)- Hp a i r -\- H res ~\~ H ro i (6) 

where H sp is a spherical one-body Hamiltonian (a spherical Nilsson model), H pair represents 
the pairing residual interaction, H res stands for the consistent with the one-body Hamiltonian 
the long-range residual interaction and H rot is the Hamiltonian of the rotor H. We use the 
ansatz of separable multipole-multipole and spin-multipole-spin-multipole forces including 
the isoscalar and isovector components 

Hres = -jEE «a[t] Ml\r}M x ,[r] - W £ «u[r] S^[t]S^[t] (7) 
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Here k\[t], ki\[t] are the multipole strength constants and 
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are isoscalar (r = 0) and isovector (r = 1) multipole fields, which are obtained from proton 
and neutron multipole operators of the form: 



M \J = (Vi a i\ R \( r ) Y ^(r)\q2V2)al ltTl a q2(T2 , 

S ill= (qi°i\ R i\( r )W® Y i(r)]\,,\(l2CT2)al 1(Ti a q2a2 , (9) 
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For symmetry reasons it is convenient to construct multipole fields of a good signature 

R^Flr}^ = rF[r] (10) 

where r = ±1 and R\ = e"'"' 7 *. The details about the symmetry properties of the operators 
can be found in |4| . After performing the signature transformation, the Hamiltonian Eq. (||) 
can be expressed through the quasiparticle creation (annihilation) a qa {a qa ) operators 

H =< HFB\H\HFB > + ]T s q £ a^v + H pair + # res + H rot (11) 

fc ci- 
in the Hartree-Fock-Bogoliubov approximation. Here e q is the quasiparticle energy. The 
pairing interaction H pair and the long-range residual interaction H res expressed in terms of 
the quasiparticle operators consist of terms proportional to combinations of type aWaW, 
craWa, «Waa and h.c. At the description of even-even nuclei the terms proportional 
to o^o^o^a and h.c. are not considered, since they create (annihilate) the odd number of 
particles. The terms proportional to cvcvcvcv and h.c are usually neglected in the descrip- 
tion of the vibrational states in the harmonic approximation (see e.g. fHU). The remaining 
terms proportional to ctWaa are treated in the random phase approximation (RPA). After 
solution of the RPA equation of motion for collective modes || the Hamiltonian Eq.flTTD can 
be expressed in terms of the RPA creation and annihilation phonons QJ V , Q v 



H RPA = £ u v {QlQ v + \) + \y. K> (12) 

and the corresponding phonon energies u u . The last term in Eq.([l^) represents the contri- 
bution of the spurious modes, for example, related to the rotation and translation in the 
coordinate space or to the rotation in the isospin space The RPA phonon defined by 

the collective coordinate X v and momentum V v has the following two-quasiparticle structure 

Ql = ^ - = E K>L-«L~ - (1 3 ) 
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Notice that Eq.(|H|) is valid only for the positive signature phonons which consist of the two- 
quasiparticle states A similar expression can be written for the negative signature 
phonons built from the two-quasiparticle states at k iOt k i and aL-a4f- Since in the phonon+rotor 
model the rotation is treated adiabatically, there is a degeneracy with respect to the quantum 
number of signature and all physical properties can be described within the space of the 
positive signature states only . In this approximation the wave function of the system can 
be written in the following form 



\vIMK") = J 16 J^ + 1 {d'mMI^k) + (-IYv^MR^k)} (14) 

where I^Wk") is the intrinsic wave function determined by one-phonon or multi-phonon 
states. The solution of the RPA equation of motion leads to the system of algebraic equations 
which defines the characteristic equation |3]14[]. The size of this equation is determined by 
the number of the operators involved in the residual interaction. It is clear that the phonon 
energies being the roots of this equation and the structure of the two-quasiparticle amplitudes 
ip u and ip u are dependent on the type of the residual interaction. In the other words, the 
presence or absence of the coupling between the electric and magnetic modes should influence 
on the structure of the RPA solutions. 
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III. DISCUSSION 



In order to avoid the cumbersome procedure of finding all phonon energies uj u needed for 
the determination of the strength function S(XX, J5 7 ), Eq.(ffl), we developed the averaging 
technique j3J. The advantage of this procedure consist of the avoiding: i)the resolving 
the many-dimensional characteristic equation for each single root ui v and ii) calculations of 
corresponding eigenvectors of the Hamiltonian and the B(XX, u u ). We consider the averaged 
strength function 



(15) 



with the averaging function such 
A 1 
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dx pa(x) 
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In Ref. @ it was shown that the results are fairly independent of the choice of the averaging 
function. Considering the function Eq . (|15"|) as a complex function with a complex argument, 
we apply the Cauchy theorem and finally obtain the following result 
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Here D is the characteristic determinant of the RPA algebraic equations and p x ~ IK is 



a quasiparticle matrix element of the transition operator. The matrix B is similar to the 
matrix D, however, it has the dimension by 1 greater than that of the matrix D 



(18) 



Sn(w) = , B a (uj) = -B u (oj) , (i = 2, . . . ,n) 
B<j(w) = A-i,i-i(^) , (*, j = 2, . . . ,n) 
The first column and the first row of the matrix B contain sums of the product of quasiparti- 
cle matrix elements of the transition operator and of the operators of the residual interaction 
(see details in [Q). The first term in Eq.([H]) is related to the residual interaction and the 
second term is determined by the mean field. 



6 



As the example of the coupling of different multipoles and spin-multipoles in the Hamil- 
tonian we consider the M2 strength function Sa(M2,E) = S«Sa(M2, E;gr — > 2~K) for 

K 

15A Sm. In Fig.l the solid line corresponds to the coupling of M\k (electric dipole) and Suk 
(magnetic quadrupole) operators. The dashed line represents the contribution of the mag- 
netic quadrupole term alone. The contribution of the plain quasiparticle term (the second 
term in Eq.(|l7|)) is shown with the dotted line. Whereas the residual interaction changes 
the distribution of the strength of the magnetic quadrupole transitions, the influence of the 
electric dipole mode to the total strength function is relatively small. It seems that in the 
region of giant resonances the coupling between electric and magnetic operators for the M2 
mode is essentially not important. The analysis of this coupling for low- lying states and for 
the other operators is in progress. 
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Figure Capture 

Fig.l The total M2 strength function vs the excitation energy E. 
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